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Tidal torques acting on a gaseous accretion disk around a binary black hole can create a gap in the 
disk near the orbital radius. At late times, when the binary inspiral timescale due to gravitational 
wave emission becomes shorter than the viscous timescale in the disk, the binary decouples from the 
disk and eventually merges. Prior to decoupling the balance between tidal and viscous torques drives 
the disk to a quasistationary equilibrium state, perturbed slightly by small amplitude, spiral density 
waves emanating from the edges of the gap. We consider a black hole binary with a companion of 
smaller mass and construct a simple Newtonian model for a geometrically thin, Keplerian disk in the 
orbital plane of the binary. We solve the disk evolution equations in steady state to determine the 
quasistationary, (orbit-averaged) surface density profile prior to decoupling. We use our solution, 
which is analytic up to simple quadratures, to compute the electromagnetic flux and approximate 
radiation spectrum during this epoch. A single nondimensional parameter g, equal to the ratio of 
the tidal to viscous torque at the orbital radius, determines the disk structure, including the surface 
density profile, the extent of the gap, the existence of an inner disk, and the accretion rate. The 
solution reduces to the Shakura-Sunyaev profile for a stationary accretion disk around a single black 
hole in the limit of small g. Our solution may be useful for choosing physical parameters and setting 
up quasistationary disk initial data for detailed numerical simulations that begin prior to decoupling 
and track the subsequent evolution of a black hole binary-disk system. 

PACS numbers; 98.62.Mw, 98.62.Qz 



I. INTRODUCTION 

Binary black hole (BHBH) mergers typically occur in 
regions immersed in gas, and the capture and accretion 
of the gas by the binary may result in appreciable elec- 
tromagnetic radiation. Following the detection of grav- 
itational waves from a BHBH merger, electromagnetic 
"afterglow" radiation could provide confirmation of the 
coalescence [ll-0|- Such electromagnetic radiation can 
also serve as a useful probe of the gas in galaxy cores 
or in other regions where mergers take place, as well as 
a diagnostic of the physics of black hole accretion. The 
timescale during which detectable "afterglow" radiation 
achieves its maximum value occurs when the gas is driven 
close to the remnant and ranges from several years to tens 
of decades in the case of supermassive BHBH systems 
with total masses of 10^ — 10^ Af©. Together with detect- 
ing the gravitational waves, observing this electromag- 
netic radiation may even provide a means of witnessing 
the birth of a quasar Q . 

There is also the possibility of detecting electromag- 
netic "precursor" radiation prior to the merger and be- 
fore the maximum gravitational wave emission [lo| . 
If the distant gas is nearly homogeneous and either at 
rest with respect to the binary ("binary Bondi" accre- 
tion) or moving ("binary Bondi-Hoyle-Lyttleton" accre- 
tion) and optically thin, the luminosity will peak at the 
end of the binary inspiral phase immediately prior to the 
final plunge [ll|. At this stage shock heating of the gas 
and turbulent magnetic field amplification is strongest. 
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The peak luminosity lasts for 5t ^ Mg hours prior to 
merger and then plummets sharply following the coales- 
cence. Here Afg is the binary mass in units of IQ^Mq. If, 
instead, the accretion takes place via a geometrically- 
thin, optically-thick Keplerian disk around the binary 
("binary Shakura-Sunyaev" accretion), there may be a 
late-time precursor brightening from tidal and viscous (or 
turbulent magnetic) dissipation in the inner disk. This 
radiation peaks on a timescale 5t ^ O.IMr days prior to 
merger and it remains high afterwards [10|. 

In this paper we focus on geometrically thin disks prior 
to disk-binary decoupling and well before any late-time 
brightening of the "precursor" electromagnetic radiation. 
Our calculations are based on a simplified, Newtonian 
prescription for a Keplerian disk in the orbital plane of a 
binary BHBH system with a low mass ratio. Similar disk 
equations have been integrated previously in time to fol- 
low the (secular) evolution of such a BHBH-disk system 
for selected cases (see, e.g., [1, Here we adopt the 

equations to solve for steady state. We then apply our so- 
lution to determine the orbit-averaged disk structure and 
electromagnetic radiation during the inspiral epoch prior 
to decoupling, when our quasistationary approximation 
is appropriate. 

Our simple treatment determines the quasistationary, 
orbit-averaged, surface density profile S(r) of the cir- 
cumbinary disk prior to decoupling, as well as the accre- 
tion rate, luminosity and approximate spectrum of the 
electromagnetic radiation. Our semi-analytic analysis, 
(i.e. analytic up to simple quadratures), serves to iden- 
tify some of the key physical parameters that determine 
these quantities. Several nondimensional parameters fix 
the overall shape of the density profile. Among these pa- 
rameters are the binary mass ratio q, assumed small, the 
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ratio of the disk scale height to radius, /i/r, also small, in 
the ring-like gap in the disk at the orbital radius of the 
secondary, and the power-law n defining the variation of 
the disk viscosity with radius, v (x r". Other param- 
eters, such as the total binary mass M, determine the 
physical scale of the disk and, together with the accre- 
tion rate M, the characteristic luminosity and frequency 
of the emitted radiation. 

We find that many of the properties of the disk, such 
as the surface density profile, the extent of the gap at 
the orbit of the secondary, the existence of an inner disk 
and the accretion rate, are determined by a single nondi- 
mensional parameter g. This parameter is essentially the 
ratio of the tidal to viscous torque at the orbital radius. 

A reliable description of the accretion flow and asso- 
ciated radiation from a merging BHBH binary really re- 
quires a radiation magnetohydrodynamics simulation in 
full general relativity in a 3 -|- 1-dimensional, dynamical 
spacetime. Such simulations have yet to be performed 
in any detail. However, Newtonian hydrodynamic sim- 
ulations incorporating some of the relevant physics have 
been performed at various levels of approximation (see, 
e.g., 0,1, andj general relativistic simulations are 
underway (e.g., [TTI. IT3l - l20l |). The model discussed here, 
although based on a simplified description, can help se- 
lect input parameters and identify scaling behavior for 
such simulations. In addition, the resulting quasistation- 
ary description of the disk provides approximate initial 
data for numerical simulations that begin prior to binary- 
disk decoupling. 

We adopt geometrized units and set G = 1 = c below. 



II. BASIC MODEL 

A. Timescales and Overview 

There are several characteristic timescales that deter- 
mine the structure of a gaseous disk in a BHBH system. 
The orbital timescale in a Keplerian disk with orbital 
angular frequency il/f (r) — {M/r^Y^^ is given by 



tor-hir) = 27r ( — 



1/2 



(1) 



where r is local disk radius. The viscous timescale is 
given by 

t ■ {r)= ^ « Mdiskr'»g(r) _ 2 r2 
(dE/dt)vis T^is 3 

where v is the shear viscosity, Tvis is the viscous torque 
and Mdisk ~ r^E(r) is the disk mass. The timescale 
associated with the gravitational tidal torque on the disk 
due to the companion is 



(dS/dt)tid 



Mdiskr-^»g(r) 



(3) 



where Ttid is the gravitational tidal torque. The inspiral 
timescale of the binary due to the emission of gravita- 
tional waves is 



tew (a) 



{da/dt)Gw 16 M^C 
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where a is the binary separation and ^ = 4(7/(1 + q)'^. 
Finally, the binary orbit is also changed by the torque 
exerted by the disk on the secondary of mass m = qM. 
The timescale is 
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The inequality torb(?') ^ ^vis('') is satisfied through- 
out the disk at all times. For most of its history, the 
disk evolves on a slow, secular, viscous timescale and 
not on a rapid dynamical timescale 2l[; viscosity main- 



tains the gas in nearly Keplerian circular orbits. In a 
quasistationary state, where the eff'ects of viscosity and 
tidal torques balance each other, we have tf^^ ~ ^vis- 
Consider the typical case where the disk mass satisfies 
Mdisk < TO. It follows from Eq. © that 4d''(«) < Cd- 
Hence, as long as the disk is quasistationary, we have 
ivis(a) ^ i™d- However, i™^ can be longer or shorter 
than tew during the early inspiral history. Prior to de- 
coupling, even as the binary separation shrinks due to 
the combined eff'ects of gravitational radiation and tidal 
torque back-reaction, the orbital separation always re- 
mains momentarily "frozen" while the disk adjusts to 
tidal-viscous torque balance to maintain quasiequilib- 
rium. Eventually, as the binary inspirals further, the 
timescale due to gravitational wave emission becomes 
shorter than the viscous timescale and the binary de- 
couples from the disk P, [^, (2^ . We note that as long as 
Afdisk/TO <C 1, ivis ^ t^d holds. Hence at decoupling and 
thereafter, the orbital decay is driven by gravitational ra- 
diation, not the torque of the disk. Quasiequilibrium no 
longer applies after decoupling. 

The orbital radius at which decoupling begins can be 
estimated from the relation tGw(a) ~ /?ivis(2Aa). Here 
A = r/2a ^ 1 is the nearly constant ratio between the 
disk edge and the orbital separation prior to decoupling 
and is determined by a balance between viscous stresses 
in the disk and gravitational tidal torques from the bi- 
nary [13, HI]- The parameter /3 ~ 0.1 roughly accounts 
for the shortening of the viscous timescale at the edge 
where the surface density E is very steep [l^j ■ An a-disk 
with a viscosity law iy{r) — (2/3)Q;Pgas/(pf^K), where p 
is the gas density and Pgas is the gas pressure, yields a 
decoupling radius ad given by [l|, 0] 

^ « 126a:;'/'°5-49/200;^7/10^2/25(^^^^^17/40^-ir/200 

. (6) 

where a ~ 0.1q;_i, (3 = 0.1/3_i, S = 37rE(ad)j^(ad)/A'/Edd 
and e = 0.26lo.2- Here MEdd = 'iirMmp/ {qar) is the Ed- 
dington accretion rate, gt is the Thomson cross-section 
for electron scattering, rj is the radiative efficiency and 
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< 1 is a porosity correction factor applied to the 
scattering-dominated optical depth [25|]. The radiative 
efficiency 77 = LEdd/M^dd may be set to 1/12 for a New- 
tonian thin disk disk with an innermost stable circular 
orbit (ISCO) at risco = 6M, the value appropriate for a 
nonspinning black hole. 

Most earlier work has focused on the late binary epoch 
following decoupling when a < ad- During this epoch 
the outer disk remains almost frozen until after binary 
inspiral and coalescence, at which point the disk diffuses 
inward and fills up the hollow (in the case of equal-mass 
black holes), or gap (in the case of a low- mass compan- 
ion), on a (slow) viscous timescale (see, e.g., d, 0| and 
references therein). Eventually the gas reaches the ISCO 
of the remnant black hole, where all torques and sur- 
face densities vanish in a first approximation. Ultimately 
the disk settles into steady-state equilibrium around the 
black hole remnant. In this paper we shall be concerned 
primarily in the early binary epoch prior to decoupling 
when a > Qd- As we have described, during this epoch it 
is the binary orbit that remains nearly frozen while the 
disk adjusts to the combined tidal and viscous torques, 
the effects of which balance each other in steady state. 
Thus, up until the BHBH orbit reaches decoupling, the 
disk evolves quasistatically as the orbit shrinks. 

Because torb('") ^ ivis(^) over the entire history of 
a quasistationary disk, the disk radial velocity satisfies 
Vr <C rflK{r) and the full hydrodynamic equations re- 
duce to secular (conservation) equations to describe the 
orbit-averaged evolution of a Keplerian thin disk [26| . 
While obtaining the spiral density wave perturbations in- 
duced by the tidal torques near the edge of the disk does 
require the full set of hydrodynamical equations, deduc- 
ing the underlying orbit-averaged disk profile does not. 
For the epoch prior to decoupling these secular equations 
can be solved in steady state to give the quasisteady disk 
structure for each value of the separation a > ad- We per- 
form this calculation below in the low-mass limit where 
q<l- 



B. Key Equations 

1. Disk Evolution 

The evolution of a geometrically thin, nearly Keplerian 
disk is determined by combining the equation of mass 
conservation. 



ld{rJ:vr) 
dt r dr 
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with the equation of angular momentum conservation 
9(Sr2f7) ld{r^Vrr^n) 1 dG 



dt 



dr 



2'Kr dr 
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to obtain an evolution equation for the surface density, 
S(i,r), 



9S _ Id 
dt 2'Kr dr 



dij^y^ dG 
dr I dr 
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the 



Here G = — Tvis + Itid is the total torque, Tvis is unc 
viscous torque, Ttid is the tidal torque on the disk from 
the presence of the secondary, and i7 = is the orbital 
frequency. The viscous torque density is given by the 
standard equation [1, [13, 



dr dr \ dr 



(10) 



We approximate the (orbit-averaged) tidal torque den- 
sity by using the expression adopted by Armitage and 
Natarajan [9| 



dTf, 



tid 



dr 
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where A(r, a) is given by 



A: 



- {fq'M/2r) [r/Apf , 
+ {fq'M/2r) {a/ Apt , 



r < a 
r > a 



(11) 
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In Eq. (|12|) / is a dimensionless normalization fac- 
tor and Ap is given by Ap = max(|r — a\,h)- Cali- 
brating the above expression for the tidal field against 
high-resolution, hydrodynamical simulations in two- 
dimensions for a low-mass, black hole secondary interact- 
ing with an outer accretion disk, Armitage and Natara- 
jan find that the value / sa 0.01 best fits the simulation 
results. Equations (fTTj) and (|12|) furnish a reasonable an- 
alytic approximation to the results obtained from sum- 
ming over the pointlikc contributions from the Lindblad 
resonances in the disk [28„i29,]. (Similar, but slightly dif- 
ferent, forms for the tidal torque also have been used in 
the literature; see, e.g. . flOl p3 C)l-l32|. For an analysis in 
general relativity, see [3j, IsJ.) Assembling the above 
expressions then yields the final evolution equation 



dt 



ld_ 

r dr 



dr 



2AS]r3/2 
Ml/2 



(13) 



The rate at which the secondary black hole migrates is 
determined both by back-reaction to the tidal torquing 
of the disk and by gravitational wave emission. 



a — atid + OGW , 



where 



4^1/2 rr. 



Otid 



(14) 



(15) 



and where oqw is given by Eq. (jH). In Eq. the 
integration is over the entire disk, although most of the 
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contribution from tidal torques arises close to the gap 
boundaries near r « a. 

The accretion rate onto the primary may be calculated 
from 



M{t,r) = 27rrS(-u, 



Or 



Combining Eqs. © and yields 



— - — — Af 
dt 2-Kr dr 



dG 
dr 



(16) 



(17) 



The coupled evolution Eqs. ([T^ and (ITil) have been 
integrated in time previously to explore select cases (see 
e.g. Here we want to solve them for quasistation- 

ary flow in general. We will thereby obtain initial data 
for the construction of disks around low-mass binaries 
prior to decoupling that is valid in general cases. Thus 
we set a = constant and 91] /dt = in Eq. (IT3)) to obtain 
the density profile. According to Eq. ([T7|) the resulting 
quasisteady accretion rate Af is independent of r and can 
be obtained by solving Eq. (I16|) once the density profile 
has been determined. 

In steady state, Eq. becomes a second-order el- 
liptic equation in r, for which we impose the following 
boundary conditions: 



b.c.'s 



(i^i;)out, r 
0, r 



^out 



(18) 



In Eq. Tout is the outer radius of the disk. We set the 
ISCO in the disk equal to risco = 6M, the value appro- 
priate for a nonspinning black hole remnant. Typically, 
''out ^ ''isco and in some cases we shall take rout — >■ oo. 

As mentioned above, the accretion rate M in steady 
state is independent of r. In steady state Eq. ([13]) admits 
a first integral which, when combined with Eq. (|16p yields 
the first-order ODE 



M = 27r 



d(ri/2i,s) 2AEr3/2 



= constant 



dr Ml/2 

(19) 

While we could solve the second-order Eq. (flSl) in steady 
state for the density profile, followed by Eq. for the 
accretion rate, it is simpler to integrate the first-order 
Eq. (dH). We do this in Sec. ImH below. 



2. Electromagnetic Radiation 

The local radiated emission from the disk arises both 
from viscous and tidal dissipation. The rate of viscous 
dissipation per unit surface area is (26j 



D.is{t,r) 



9 M 



(20) 



The integrated rate of tidal dissipation may be obtained 
from the change in the binary binding energy due to tidal 
forces. 



4id(t) = 



2a 



(21) 



tid 



where (a/a)tid is given by Eq. We follow [l^ and 

assume that tidal dissipation results from the local damp- 
ing of spiral density waves that mediate the binary-disk 
interaction and is thus proportional to the local genera- 
tion of these waves, 



r, u \ ^ TP \dTtid/dr\ 

47rr dr\dTtid/dr\ 



(22) 



Approximating the emission as thermal blackbody ra- 
diation, the local disk surface temperature Ts(i, r) may 
be equated to the effective temperature, which is deter- 
mined from the total dissipation rate D = Dyis + ^tid 
according to 



aTfit,r)=D{t,r) 



(23) 



where a is the Stefan-Boltzmann constant [35j. 

Given the surface temperature, the quasistationary 
specific flux F^, (t) measured by an observer at distance d 
whose line of sight makes an angle i to the normal to the 
disk plane is determined by integrating over the entire 
disk surface, 



2tt cos i 

d2 



BjTs{t',r))rdr, (24) 



where B^{Ts{t' ,r)) is the Planck function, t' = t — d is 
retarded time and v is the photon frequency [36| . Equa- 
tion (|24|) is best evaluated in terms of a nondimensional 
function f*(t',x) of nondimensional frequency x accord- 
ing to \6] 



27r cos n5 aT^ 

d2 TT^ 



^isco/ ' •^)' 



where 



du u- 



exp (xT^/Ts) — 1 ' 

and where we have introduced the parameters 

aT^ = 3M Mrcm/87rr? hi^, = kT, , 
X = hv /kTf, = V /v^,, u = r/risco- 



(25) 



(26) 



(27) 



Here Aficm = 37r(z^I])out is the stationary accretion rate 
of an infinite disk onto the remnant black hole follow- 
ing merger [see Eq. ([3T|) ]. The quantity T* provides a 
convenient estimate of the characteristic temperature in 
the main radiating region near rjsco of the final equilib- 
rium disk following merger, and hu^, is the characteristic 
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frequency of the emitted thermal radiation from this re- 
gion. The specific luminosity L^(t), summing over both 
surfaces of the disk, is related to {t) according to 



cos I 



(28) 



The total luminosity L{t) integrated over all frequencies 
is then given by 

/"OO pf\ fOO 

m= dvL,{t) = —aTtrl^J f* {f ,x)dx {2^) 
Jo Jo 

Given the quasistationary density profile, we wish to 
perform these quadratures for the quasistationary elec- 
tromagnetic spectrum prior to decoupling. Equation (j29p 
yields 



Z^visSTrrdr + Et 



id 



(30) 



We note that the final stationary equilibrium disk fol- 
lowing merger yields well-known analytic density and 
temperature profiles, as well as analytic integrated fluxes 
and luminosities (see, e.g., @), and these quantities pro- 
vide useful checks on the numerical quadratures. For 
example, after the merger, the tidal torque is absent and 
the equilibrium density profile is given by the familiar 
Shakura-Sunyaev result ([131; see also [la, [13, HI] and 
references therein) for a thin disk around a single black 
hole remnant. 
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The corresponding luminosity is given by 
MM 



(31) 
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^ 3/2 



2n, 



[post — merger] 

(32) 



where 



and where 



2fq'M^/Vj^l 



t'out 



gsi s > si 
—gs^ s < si 



(35) 



(36) 



Evaluating Eq. ([34]) in steady state yields the second- 
order ODE 



ds 



ds 



max(]s2 — s1\,s'^h) 



2/ =0 



(37) 

which must be solved for s G [s2, 1] subject the boundary 
conditions 



b.c.'s : y = y = 



1, s = 1 

0, S — S2 



(38) 



Equation ((37|) admits a first integral, just as in Eq. (|16l) . 
so that we need only solve the first-order ODE 



ds 



where 



5lf) 



max(]s2 — sfl.s'^h) 



(39) 



(40) 



III. STRUCTURE EQUATION: STEADY-STATE 
SOLUTION 



A. General Solution 



To solve Eq. p9l) we introduce the function 



F{s) = / /(s')ds' 



(41) 



3. Structure Equation: Nondimensional Units and Scaling 



To solve Eq. (1131) in steady state for general cases, and 
to help identify scaling behavior, it is convenient to in- 
troduce the following nondimensional variables: 

s = (r/rout)^^^, si = (a/rout)^^ J S2 = (nsco/»'out)^^^, 
E = E/Eout, v = i'/yout, y = sT,, y = h'y, 
h = h/r, T = i/2ivis(?'out), m = M/ (37ri'outSout)- (33) 

In terms of these variables, equation (fT3)) becomes 



dy _ 1 d-'jPy) 1 d J 



1 



max(]s2 — sl\,s'^h) 



(34) 



The ODE can be rewritten as 
d 



ds 



y) 



and is readily integrated to give 



y{s) = e-^(^) 



b-rh e^^^'^ds' 



(42) 



(43) 



The boundary condition y(l) — 1 gives 6=1, while the 
condition y{s2) = gives 



m 



(44) 
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Hence the solution is 



vis) 



(45) 



To compute rfi and y{s) numerically, it is convenient 
to introduce the function 



Gis) 



Functions F and G satisfy the coupled ODEs 

F'is) = -f{s), 

G'{s) = -e^^'\ 

with the initial conditions F{1) — G(l) = 
tions (|44p and (pSj) can be written as 



G{S2) ' 



1 



Gis) 



Gis2) 



(46) 

(47) 
(48) 
0. Equa- 

(49) 
(50) 



In some cases, G{s) may become large when s ap- 
proaches si. To treat this complication we introduce the 
function H{s) — ln[G(s)], which satisfies the ODE 



H\s) 



G{s) 



^Fis)~H{s) 



(51) 



In practice, the coupled ODEs (j47|) and (|48l) are inte- 
grated from s = 1 to s = s*, where si < s* < 1. Then 
Eqs. (|47l) and ([5T|) are integrated from s ~ s* to s ~ S2 
with the initial condition H{s*) = ln[G'(s*)]. Having 
computed the functions F and H for s € [s2, 1], irt and 
y{s) are obtained by 



vis) 



-H{S2) 



Dis) 



1 



^H(s)-H(s2) 



(52) 

(53) 



We shall study the solution by numerically integrating 
the above equations for different choices of parameters in 
Section HireH below. 



B. Limiting Cases and Asymptotic Behavior 

Before obtaining numerical solutions for general cases 
it is instructive to evaluate Eq. (15^ analytically for the 
disk structure in limiting regimes. We first observe that 
the parameter g defined in Eq. ([55|) may be evaluated as 



9 = 2nfq^ 



torh / out 



(54) 



where we have used Eqs. ([T]) and ([2]). We also find it 
useful to introduce another nondimensional parameter 



2h^ Uorb 



■iv(a)(h/af ■ 



(55) 



The physical meaning of g can be understood as fol- 
lows. Consider a ring of disk material of width h near 
r — a. The tidal torque action on the ring is [see 
Eq. PT|) ] Ttid(a) ~ 27rASa/i. The viscous torque is 
rvis(a) = -2Tra^i^{a)j:dMK{a) [see Eq. Hence 



rtid(a) _ 2h\A{a)\ 
rvis(a) 3i^(a) 



1/2 



9- 



(56) 



The parameter g therefore measure the relative impor- 
tance of the tidal to viscous torques at r = a: when g is 
small the tidal torque is unimportant, otherwise it plays 
a significant role. The variables /, q and h in Eq. ((55)) 
are each always less than unity, while the ratio ivis/^orb is 
always larger than unity; hence g can vary between zero 
and infinity, depending on the model (the same conclu- 
sion also applies to the parameter 17). 



1. Negligible Tidal Torques 

In the limiting regime g ^ tidal torques are negligi- 
ble and Eq. (|39p. may be solved together with boundary 
conditions ([55]) to yield 



y ^ m{s - S2), (ff = 0), 



and 



m 



I-S2 



Restoring units, Eqs. (|57)) and ([58)1 translate to 



M r 
3^ 



,1/2 /^l/2 
isco/ 



(57) 



(58) 



(59) 



and 



M 
37r(iyE) 



1 



1 _ ^1/2 /^l/2 
isco/ out 

1, r-out 00, 



(60) 



respectively. As expected, when tidal torques arising 
from the companion are negligible, Eqs. (|59p and ([50)) 
reduce to Eq. ((3T|l , the quasistationary profile and accre- 
tion rate for a disk around a single black hole. 



2. Strong Tidal Torques 

In another limiting regime g 3> 1, tidal torques are 
strong and we will see in the following analysis that the 
accretion is shut off. 

In the thin-disk limit /i — >■ 0, the function /(s) changes 
from a large negative value to a large positive value at 
s ~ si. The function F{s) defined in Eq. (|¥T|) thus has 
a narrow peak near s = Si. Hence we may approximate 
TO in Eq. (|44p by integrating over the small region 
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around s = Si where |/(s)| reaches maximum. We first 
approximate f{s) near s = Si by a step function 




Si < s < s~^, 

< S < Sl, 



(61) 



where s' 



1 



Fis) 
where 



- si/\/l - h « + h/2) and = 

si(l - h/2). Using Eq. gl]) we find 



2(s-si) 



Sih 



for 



< s < 



(62) 



fis)ds 



P(si) 



(is 



2\4 



= |[l + 0(/.)]. (63) 
The accretion rate m can then be approximated by 



1 

m 











/ (is exp 


.9 - 2 


s — Sl 


)] 






Sl/l 





,9/3 



£1^,9/3 



3 

For ^ 1, we have 



1) 



9_ c-4g/3. 
Sift. 



(64) 



(65) 



Hence when the tidal torques are strong, the accretion 
is shut off. Our numerical calculation confirms that 
Eq. is accurate when ^ 1 and /i <C 1. 

In Sec. lIIIBTl we see that m « 1 for weak tidal torques 
{g ^ 1). The above analysis shows that m « for 
strong tidal torques (e^ ^ 1). The critical transition 
between these two regimes occurs at g = ^c, which can 
be estimated by the condition 



_5^g-4Sc/3 
Sift 



1. 



(66) 



For a = lOOM, rout = lO^M (si = 10'^/^) and h = 0.1, 
Eq. ([5^ gives (jc 6. We find in Sec. lIII B"4l below that gc 
lies between 5 and 10 for a power law viscosity i/(r) cx r", 
which is consistent with this simple analysis. 

The density profile in the strong tidal torque regime is 
very different from that in the weak tidal torque regime. 
We will see in Sec. IIIIB 31 below that in the asymptotic 
region in which r 3> a, the density profile is given by 



S(r) 



-1/2 



\rout 



-(n+l/2) 



(67) 



for a power law viscosity. To estimate E near r = a, we 
first note from Eq. that G(s) < G{s2) for s > S2 
(r > nsco)- Hence y{s) < e^^^'^'> /D{s). For s" < s < 
s+, Eqs. ^ and ^ give 



S(r) 



< 



1 



sv{s) 



exp 



'9 



2(s-si) 



Sih 



< 1 (68) 



for large g. For s < s , we combine Eqs. p6)) . (|49l) and 
(l50t to write 



(69) 



It follows from Eqs. g?]), dM]) and ([Ml) that F(s') - 
F{s) < for s' < s < Sl. Hence 



S(r) . / r 

^out \ f'out 



1 



< 1 (r < a) (70) 



for small m (large g). 

We therefore see that the density plummets for r < a 
when the tidal torques are strong. This behavior can 
be understood as follows. The tidal torque Ttid pushes 
the disk matter radially outwards for r > a and inwards 
for r < a; whereas the viscous torque Tvis always drags 
the matter inwards. The tidal torques are strongest near 
r — a, which tend to suppress the inflow and create a 
gap there. When Ttid is sufficiently strong, the inflow 
is strongly suppressed near r — a. Inside the binary's 
orbit (r < a), both Ttid and Tvis push matter inwards, so 
that when Ttid is strong matter cannot accumulate in the 
inner region. Hence in steady state, there is practically 
no inner disk and the density of the outer disk near r — a 
is also very small. 



3. Asymptotic Behavior 

As shown above, when tidal torques are present, the 
accretion rate may be much lower than the value quoted 
in Eqs. ([58]) and ((60l) for accretion onto a single black 
hole. The rate m depends sensitively on g and requires 
the full numerical solution to determine accurately. But 
in the asymptotic region r 3> (3 it is always possible to 
express the disk structure in terms of m even when tidal 
torques are present, as we now show. 

Whenever 1 > s > si, Eqs. ([39]) and (j40|) reduce to 



ds 



y 



Consider the asymptotic regime where 

^(^)'«1, 
V \ s / 



(71) 



(72) 



in which case the second term in Eq. ([7T|) is typically 
much less than the first term and therefore can be 
dropped. The solution to the resulting equation is 



y 



ms - 



1 



TO, 



(73) 
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m + (1 — m) 



-1/2 



(74) 



where we have used the boundary condition y(l) = 1 
(Eq. [35]). As we have argued in Sec. lIIIB^ whenever the 
tidal torque parameter satisfies g ^ ffc we have m ss 1, 
while whenever g ^ gc we have m ^ 1. Not surprisingly, 
whenever tidal torques are unimportant, the asymptotic 
density profile refiected by Eqs. (1751) and (fMl) is the 
same as the asymptotic profile given by Eqs. ([57| and 
(|59|) (i.e., vTi w constant (z^I])out ~ M/Stt), but has a 
very different fall-off otherwise. 

As an application, take the disk viscosity to be a power- 
law profile, 

i^{r) cx r", n^ -7/2, (75) 
or P{s) = s^" with n ^ -7/2. Then Eq. ^ becomes 



as 



where 



q{s) = exp 



{2n + 7)s 



2n-l V s y 



Consider the asymptotic region where 

, 8 

(2n-(-7)s2"-- \ s 

or 



s > 



2n-l \ l/(2n+7) 

2n + 7 



Si . 



(76) 



(77) 



(78) 



(79) 



For example, if n — 1/2, our fiducial case below, Eq. ([79 
reduces to 



s > 



(!) 



1/8 



Si . 



(80) 



In the asymptotic region defined by Eq. ()79p . q{s) « 1 
and Eq. ([75]) again can be integrated to give Eq. ([75|. 
which now yields 



„2n+l 



rout 



+ (1-to) 



Tout 



-(n+1/2) 



(81) 



Numerical Integrations 

Here we adopt the numerical recipe outlined in Sec- 
tion IIII Al to evaluate the quasistationary density pro- 
file for different disk parameters. We adopt a power law 
viscosity i/(r) oc r". We note that in their vertically- 
integrated hydrodynamical simulations ref. [12] adopts 
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FIG. 1. Accretion rate m as a function of the tidal torque 
parameter g for various values of the viscosity power-law index 
— 1 < n < 2. Results are shown for binary orbital radius 
a = lOOA/, disk thickness h/r = 0.1 in the vicinity of the disk 
edge at r ~ a (see Eq. I12p . disk outer radius rout = 10^ M 
and disk ISCO at risco ~ 6M. Solid (black) line represents 
the case for n = — 1, dotted (red) line for n — 0, dashed (blue) 
line for n = 1/2, long-dashed (magenta) line for n — 1, dot- 
dashed (green) line for n = 3/2 and dot-long-dashed (cyan) 
line for n = 2. 



n = 1/2 for the value of the viscosity power-law index, 
while in their integrations of the radial secular evolution 
equations ref. [1] adopts n — 3/2, while ref. Q considers 
n = 0.4 and n = 0.065. In principle, the viscosity depen- 
dence, as well as the disk thickness, can be determined 
self-consistently once a viscosity law is adopted (e.g. an 
a-disk prescription). But due to the uncertainty in this 
law, we treat these as free parameters and postpone a 
self-consistent calculation for a future investigation. 

In Fig. [1] we show the dependence of the accretion rate 
rh — A//(37rfoutSout) on the strength of the tidal torque 
parameter g defined in Eq. (j55p for fixed binary radius 
a = lOOM, disk thickness h/r = 0.1 and disk outer- 
radius Tout = lO^Af. In Fig. [21 we show the g dependence 
of rh for n = 1/2 and n = 3/2 and for various values of 
disk thickness h/r between 0.01 and 0.15, keeping other 
parameters fixed. 

The key result emerging from these two figures is that 
for small tidal torques 5 ^ ffc the accretion rate is very 
close to the rate onto an isolated black hole for the same 
asymptotic disk parameters, while for strong torques 
g > jjc, it is sharply reduced. This behavior is predicted 
from the simple analysis in Sees. IIII B ll and rill B 21 Since 
the parameter g varies as g^, the accretion rate thus de- 
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100 



FIG. 2. Accretion rate m as a function of the tidal torque 
parameter g for various values of disk thickness h/r. Results 
are shown for the viscosity power- law index n — 1/2 (top) 
and n = 3/2 {bottom), with binary orbital radius a — lOOM, 
disk outer radius rout ~ lO^M and disk ISCO at risco = 6M. 
Solid (black) line represents the case for h/r = 0.01, dotted 
(red) line for h/r = 0.05, dashed (blue) line for h/r — 0.1, 
and long-dashed (magenta) line for h/r = 0.15. 



creases as the binary mass ratio increases. If we define 
(jc to be the value of g so that rh = 0.5, we find gc lies 
between 5 and 10 for all the cases considered here, which 
is consistent with the value 6 estimated in Sec. IIII B 21 

In Fig. 13] we show the surface density profile for dif- 
ferent tidal torque strengths g, fixing the viscosity index 
to n — 1/2, disk thickness h/r — 0.1, binary separation 
a = lOOAf and rout = lO^M. It is evident from the figure 
that when the tidal torque is large and g ^ 5, an appre- 
ciable gap develops in the density profile at the orbital 
radius of the secondary near r = a. Moreover, the gap 
widens and the density inside the gap falls off sharply 
as g 15. By contrast, when the tidal torque is small 
and ^ ^ 5, the gap shrinks and the density profile ap- 
proaches the stationary profile around a single black hole 
of the same mass as the binary (Eq. ICTl) . We thus find 
explicitly that the stationary solution yields an inner as 
well as an outer disk for sufficiently weak tidal torques 
(e.g., binaries with small q) and hardly any inner disk for 
strong tidal torques (e.g., binaries with large q), as pre- 
dicted in Sec. IIII B 21 Correlating the values of m with 
the density profiles for E for each value of g shows that 
the asymptotic profiles plotted in Fig. |3] agree with those 
derived in Eq. ([5T|) . 

We note that for large g when the tidal torque is very 
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FIG. 3. The surface density profile S plotted for different 
values of the tidal torque parameter g. All results assume 
h/r = 0.1, n = 1/2, a = lOOA/, rout = lO'^M and n^co = 6M. 
Top: g = 5000 (black solid line), g = 500 (red dotted line), 
g = 50 (blue short-dash line) and ^ = 15 (magenta long-dash 
line). Middle: g = 5000 (black solid), g — 15 (magenta long- 
dash line), g — 5 (green dot-short-dash line), g — 2.5 (cyan 
dot-long-dash line), and ^ = (orange long-dash-short-dash 
line). Bottom: Same as the middle figure but plotted in linear 
scale in r/M. 
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strong, there is essentially no inner disk and the gas is 
unable cross the binary orbit. The surface density S 
plummets near r — a. Hence the disk profile can also be 
obtained by imposing the boundary condition E = at 
r = a for sufficiently large g. We have verified that for 
large imposing the boundary condition at r = a gives 
the same disk profile as that obtained by imposing the 
boundary condition at r = t'iscOj confirming our intu- 
ition. Since we do not know a priori exactly how high 
g must be for this approximation to apply, we must im- 
pose the true physical boundary condition at r = nsco 
in general. 



IV. ELECTROMAGNETIC SPECTRUM 

Here we compute the (approximate) electromagnetic 
spectra from stationary disks for some representative, as- 
trophysically plausible systems. We are particularly in- 
terested in probing disks at the time of decoupling and 
immersed in gaseous regions that will yield stationary 
post-merger accretion rates onto the remnant black hole 
near the Eddington limit. After decoupling, our station- 
ary disk solution is no longer valid. To specify a model 
we first choose the parameters M, a/M, rout A/ ^ 1 
and n. We determine the decoupling separation by 
the condition tGw(ad) = l^Uisi^ad) (setting /3 = 0.1), 
which yields 



'ad\2 ^ 128Af/3C 



128Af/3C /rout 
15j^out V 2ad 



or 



M 



128M/3C /routV 
ISt^out V2mJ 



l/in+2) 



(82) 



(83) 



We next set h = h/r = 0.1 ioi the disk thickness near r = 
a. To establish the scale for the density and disk size in 
physical units we fix Eout and rout , which determines the 
disk mass Mdisk- Finally, we determine Uo^t by specifying 
the accretion rate onto the black hole remnant, Mj-om = 
S'Ti'outSout (see Eq. [Sniin the limit risco/^out ^ 1) to be 
a fraction 7 of the Eddington value, whereby 



MEdd 



S'Tt'outSout 
MEdd 



which gives 



CTT^out 



(84) 



(85) 



The above choices determine the torque parameters g 
and g according to Eqs. (|35|) and ([55l) . setting / = 
0.01. These choices specify a unique, albeit approximate, 
model for a stationary Newtonian disk about a binary 
black hole system and determine the disk surface den- 
sity and the gas accretion rate. Given such a model, the 
electromagnetic spectrum can be computed following the 
prescription outlined in Section III B 21 above. 
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FIG. 4. Surface density of disks around a binary black hole 
just before decoupling. Shown here are cases for the binary 
mass ratios (a) g = (black solid line), (b) q = 2 x 10~^ 
(red dotted line), (c) g = 3 x 10^'^ (blue short-dash line), (d) 
q = 4: X 10~^ (magenta long-dashed line), (e) g = 5 x 10~^ 
(green dot-short-dash line), (f) q = 10~^ (cyan dot-long-dash 
line), and (g) q — 10^^ (orange long-dash-short-dash line). 



A. Numerical Integrations 

We illustrate the calculation by considering a binary 
black hole with a IO^Mq primary. We take the viscosity 



power-law index to be n = 0.5, and set ro 



— in3 



lO-'M, 



7 — 0.1, and Eout = 1-5 x 10 gcm~^. The decoupfing 
separation aj. is calculated from Eq. ([83]) . Our choice 
of parameters gives Afdisk ^ 1O'^M0 for < g < 0.1 
(see Table |T| , which is much smaller than the mass of 
the secondary BH for q > 10^^. These values for the 
disk density and mass are comparable to those consid- 
ered in . We evaluate the disk profiles and elec- 
tromagnetic structure for different choices of mass ratio 
g < 0.1. Table U summarizes our results. Scaling behav- 
ior for different choices of binary and disk parameters 
is presented in Appendix lAl We see that the decoupling 
separation oc [see Eq. (|83l) ] increases from 12. 6M 
for g = 2 X 10^3 to 56M for q = 0.1. The accretion 
rate drops rapidly towards when q > A x 10"'^, which 
corresponds to g > 6.3, and is consistent with Fig. [TJ 

Figure [4] shows the surface density of the disk for dif- 
ferent values oi q. As in Fig. [31 for small q, there is a 
dip near r = a in the density profile, but an inner disk of 
substantial density is present. When q > 4 x 10~^, the 
accretion rate drops to a small value and E falls sharply 
near r = a, leaving essentially no inner disk. We also 
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TABLE I. Properties of disks around a binary black hole at the decoupling separation. The mass of the primary black hole 
is set to AI = IO^Mq, outer disk radius rout = 10'^ Af, surface density Eout = 1-5 x 10*gcm~^ and the Eddington parameter 
7 = 0.1. 



Case 


q 




g 


torb (days) 


tew (yes) 


aa/M Mdisk(MQ) 


m 


M{Mq yr' 


-1) L(ergs-i) 


a 
















660 


1.08 


0.29 


1.3 X lO"*^ 


b 


2 X 10" 


'3 


1.6 


1.6 


16 


12.6 


660 


1.07 


0.29 


1.5 X lO"*^ 


c 


3 X 10" 


3 


3.5 


2.0 


20 


14.8 


700 


0.953 


0.26 


3.0 X lO"*^ 


d 


4 X 10" 


3 


6.3 


2.4 


23 


16.6 


890 


0.381 


0.10 


9.1 X lO"*^ 


e 


5 X 10" 


3 


9.8 


2.8 


27 


18.2 


1000 


0.0197 


5.3 X 10" 


^ 1.1 X lO''^ 


f 


10-2 




39 


4.2 


40 


23.9 


1000 








6.5 X lO"*^ 


g 


10-^ 




3900 


15 


145 


56.0 
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FIG. 5. Ratio of tidal heating to viscous heating Dtid/fvis 
as a function of radius for disks around a binary black hole 
just before decoupling for binary ratios (b) g = 2 x lO""^ 
(red dotted line), (c) g = 3 x 10"^ (blue short-dash line), (d) 
g = 4 X 10"^ (magenta long-dashed line), (e) 5 = 5 x 10~^ 
(green dot-short-dash line), (f) q — 10~^ (cyan dot-long-dash 
line), and (g) q = 10~^ (orange long-dash-short-dash line). 



FIG. 6. Electromagnetic spectrum of disks around a binary 
black hole at redshift z just before decoupling. Shown here 
are cases for the binary mass ratios (a.) q — (black solid 
line), (b) g = 2 X 10"^ (red dotted line), (c) g = 3 x 10"^ (blue 
short-dash line), (d) q — Ax 10"'^ (magenta long-dashed line), 
(e) g = 5 X 10"^ ( green dot-short-dash line), (f) g = 10 
(cyan dot-long-dash line), and (g) g = 10~^ (orange long- 
dash-short-dash line). 



see that the density in the asymptotic region changes 
from r^^ for small m (large q) to r"^/^ for substantial m 
(small q), as predicted by Eq. ([5T|) . For a fixed Sout, the 
steeper decrease of density in the small accretion cases 
(g > 4 X 10~^) give larger S in the asymptotic region, 
which results in a larger disk mass as shown in Table HI 
However, as q increases further, the decoupling radius ad 
moves to a larger radius so that the inner radius of the 
disk moves out as well. This behavior leads to the small 
decrease in disk mass for q > 10^^ as shown in Table HI 
Figure [S] shows the ratio of tidal heating to viscous 
heating i'tid/i'vis as a function of radius for each case. 



In all cases with q ^ 0, the effect of tidal heating dom- 
inates over the viscous heating at radius r a, but de- 
creases rapidly both at larger and smaller radii. This 
result is expected since the tidal torque is created by the 
gravitational force of the secondary black hole, which de- 
creases with the distance from the secondary. It can 
be proven easily from Eqs. (gl]), (HI]) and ^ 

that Dtid/Dvis ^ r^lAl/v. Hence Dtid/Dvis has a peak 
near r = a, -Dtid/-Dvis ~ r'l'^ jir — a)* for r ^ a and 

Figure [S] shows the electromagnetic spectrum of the 
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disk for each value of q. We see that the himinosity in- 
creases substantially at first as q increases above zero and 
the accretion rate m drops. This is due to the large in- 
crease in E as shown in Fig. U) The peak frequency also 
shifts to a higher value because the enhanced viscous and 
tidal dissipation give rise to higher temperatures in the 
disk. For 9 ^ 5 x 10""^, however, the tidal torque of the 
binary becomes large enough to expel the disk material 
inside its orbit and shut off the accretion. As a result, the 
edge of the disk moves to a larger radius where the tem- 
perature is lower. This effect shifts the peak frequency 
back to a lower value. Since a large fraction of the electro- 
magnetic radiation comes from the innermost part of the 
disk where the temperature is high, the disk luminosity 
drops when g ^ 5 x 10~^ (see Table as the inner edge 
of the disk moves to a larger radius. We see from Table U 
that the total luminosity of the disk for (7 = 0.1 drops to 
about the same value as the standard Shakura-Sunyaev 
disk. 

Most noteworthy is the fact that even as the gas accre- 
tion rate is shut off by the companion, tidal heating serves 
to maintain a high luminosity from the disk, comparable 
in magnitude to the luminosity from a single black hole. 
For high q cases in which there is no inner disk, the lumi- 
nosity plummets after decoupling over a timescale tew, 
as the inspiral accelerates and tidal heating drops. Later 
on the luminosity from the disk about the remnant then 
increases on a slow viscous timescale, as the outer disk 
diffuses inward [J H, 0] ■ This luminosity achieves values 
comparable to those prior to decoupling, but at higher 
frequency, since the main radiation region outside the in- 
ner edge of the disk moves inward from ad to nsco, where 
temperatures become hotter. For low q cases in which 
there is a substantial inner disk the change in luminosity 
and the frequency spectrum following decoupling is much 
less pronounced. 

Periodicities in the luminosity are expected on an or- 
bital timescale during the inspiral. This timescale is on 
the order of days at decoupling for the cases considered 
here (see Table 

Figure [B] indicates that most of the radiation emitted 
by the disk at decoupling is in visible and near infrared 
wavelengths. For a source at redshift z = 1, the appar- 
ent magnitude is about 19 and radiated in infrared wave- 
length. Assuming the disk is not obscured by interstellar 
dust, it could be observable by the James Webb Space 
Telescope (JWST) and Large Synoptic Survey Telescope 
(LSST). 



V. SUMMARY 

We have studied the effects of tidal torque on a qua- 
sistationary disk around a binary black hole with small 
binary mass ratio before decoupling. We find that the 
density profile and accretion rate is sensitive to the mass 
ratio q through the tidal torque parameter g defined by 
Eq. (1551) . For small g, the density profile and accretion 



rate are slightly modified from the standard Shakura- 
Sunyaev disk. Specifically, the density inside the binary 
orbit (r < a) is smaller than the standard Shakura- 
Sunyaev disk but is about the same in the outer region. 
A substantial drop in density is seen near r — a due to 
the strong tidal torque in that region. As g increases, 
the tidal torque becomes stronger and the density in the 
inner disk, as well as the accretion rate, decrease fur- 
ther. On the other hand, the density in the outer disk 
increases and a gap (inside of which the density is much 
smaller) develops between the inner and outer disk. The 
gap widens with increasing g. When g exceeds a critical 
value cjc, the tidal torque is strong enough to effectively 
expel the material inside and near the binary orbit and 
practically shuts off accretion. The asymptotic density 
profile also changes [see Eq. (I5T]) ]. For a power law vis- 
cosity cx r", the critical value \s gc ^ 1 and depends 
weakly on the index n and disk thickness h/r (see Figs.[T] 
and HI . 



We compute the luminosity and electromagnetic spec- 
trum emitted by the disk due to viscous heating and tidal 
dissipation at the time of decoupling. The results are 
summarized in Table U and Fig. [5] for disks surrounding 
a primary black hole of mass M = IQ^Mq, outer disk ra- 
dius r = lO'^M, surface density Eout = 1-5 x 10*gcm~^, 
and the Eddington parameter 7 — 0.1. The disk lumi- 
nosity at first increases with q. This is caused by the 
increase in density of the outer disk. The radiation also 
shifts to higher frequencies, as the enhanced viscous and 
tidal heating give rise to higher disk temperatures. The 
disk luminosity reaches a maximum near q — 5 x 10~^, 
when the accretion is shut off by the strong tidal torques. 
When q > 5 X 10^^, the decoupling binary separation 
moves to a larger radius. The strong tidal torque expels 
disk material inside and near the binary orbit, causing 
the radiation region to move to larger radii, where the 
temperatures are lower. The disk luminosity drops and 
the radiation shifts to a lower frequency. We note that 
even as the gas accretion rate is shut off by the compan- 
ion, tidal heating serves to maintain a high luminosity 
from the disk, comparable in magnitude to the luminos- 
ity from a single black hole. Most of the radiation is 
emitted in the near infrared and visible wavelength and 
may exhibit periodicities on the binary orbital timescale. 
For a source at redshift z = \, the apparent magnitude 
is about 19 (provided the disk is not obscured by the 
interstellar dust) and could be detected by JWST and 
LSST. 
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Appendix A: Scaling in disk radiation 

In this appendix, we rewrite the equations in Sec. lIIB 2l -Dtid 



Hence Eq. ([22]) becomes 

SMi^outSoutil - m{l - S2)] 



\fi-sMs) 



in terms of the nondimensional variables introduced in ^(routci) [1 2i/(si) + m(2si S2 

Sec. EHl and then establish scaling relations in terms of ^j^^ ^^^^j dissipation rate is 
the mass of the primary black hole M and the Eddmgton 
parameter 7. 



1)] 



We first rewrite Eq. ([20|) as 



D^is(t,r) = 3— 



The torque function A can be written as 



4r ' t 



(A2) 



Hence the rate of change in binary binding energy is [see 
Eq. da])] 



37rM^/routt'outSout 



tid 



a3/2 

37rAf^routfoutSout 



w 



3/2 



f{s)y{s)ds 



l-m(l-S2)]. (A3) 



The tidal torque density Eq. PT|) becomes 



dTtid _ 37rEoutJ^out f{s)y{s) 
rout s 



dr 2 
Consider the expression 



dr\dTtid/dr\ =^ 2rout / s\dTtid / dr\ds 
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(Al) where 
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3y 



1 - m(l - S2) 



1 ~ 2y{si) + m{2si - S2 - I) s^sf 



The temperature ratio is given by 
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(A9) 



3MMrom/87rr3^„ 
Gathering all the relevant formulae, the accretion rate, 
luminosity, and disk mass can be expressed in cgs units 
as follows: 



M = O.27A/0 yr" 
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7\3/4/ M ^^Z" 



(AlO) 



(All) 



i/h. = 2.91 X 10 



L = 1.26 X 10"*^^ 
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-1/4 
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D{s)ds, 



Hz, (A12) 

erg s-\ (A13) 
(A14) 
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f{s)y{s)ds 



and 



A/disk = 2.06 X lO^Afo 2 ; V 103. f y 

1.5 X lO^gcm^"^ / VIO^A// 



/ Af 



= STrEouti'out \/ A/rout x 

[l-2y(si)+m(2si-S2-l)]. (A5) 
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